Semiconductor effective charges from tight-binding theory 
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We calculate the transverse effective charges of zincblende compound semiconductors using Har- 
rison's tight-binding model to describe the electronic structure. Our results, which are essentially 
exact within the model, are found to be in much better agreement with experiment than previous 
perturbation-theory estimates. Efforts to improve the results by using more sophisticated variants 
of the tight-binding model were actually less successful. The results underline the importance of 
including quantities that are sensitive to the electronic wavefunctions, such as the effective charges, 
in the fitting of tight-binding models. 
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The Born effective charges e^, also known as transverse 
or dynamic effective charges, are the fundamental quan- 
tities which specify the leading coupling between lattice 
displacements and electrostatic fields in insulators.a In 
general the effective charges are site-dependent tensors. 
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where Pi is the polarization in cartesian direction i, and 
M^'^ is the displacement of sublattice I in cartesian di- 
rection j. However, for compound semiconductors of the 
zincblende structure, which are the focus here, it is easily 
shown that the effective charges are scalars, and are equal 
and opposite for cation and anion; it is conventional to 
use the positive cation effective charge to characterize a 
given compound. The effective charges for a variety of 
zincblende semiconductors have been computed using ab- 
initio density-functional linear-response theory ajid have 
been found to agree very well with experiment .Ha How- 
ever, it is interesting to inquire whether more approxi- 
mate schemes can give a good accounting of the effective 
charges in compound semiconductors. If so, additional 
insight into the chemical and physical factors that affect 
the might be obtained. 

One particularly attractive and well-known approx- 
imate scl^ftHie is the universal tight-binding model of 
HarrisonOtllt provides a straitforward and computation- 
ally efficient approach to calculating electronic properties 
of solids using a minimal orthogonal sp^ basis set, with 
the Hamiltonian limited to the on-site and nearest neigh- 
bor terms. The on-site elements and ep are taken from 
calculated free-atom term values, while the interatomic 
elements {Vssa-, Vspa, Vppa and VppT^) are taken to be 
species-independent "universal" constants times the in- 
verse square of the distance. Given its simplicity, the 
model is impressively successful in estimating nmoy elec- 
tronic properties of a wide variety of materials.QB 

It is thus natural to ask what the Harrison tight- 
binding model would predict for the effective charges 
of the zincblende compound semiconductors. Oddly, 



this question does not appear to have been answered 
previously. The only previous work of which we are 
aware made use of a two-center perturbation approj 
mation to obtain estimates of the effective charges 
This approach used an expedient division of the effec- 
tive charge into "static" and "transfer" charge contribu- 
tions, with theJnterpretation of the latter being open to 
some question.El The purpose of this Report is to present 
essentially exact calculations of the transverse effective 
charges computed for zincblende II- VI, III-V, and IV- 
IV semiconductor compounds using the Harrison tight- 
binding parametrization. While the results could have 
been obtained using linear-response techniques, we found 
it simpler to to compute the e^'s instead from finite dif- 
ferences, calculating directly the change in bulk polar- 
ization from a small displacement of one sublattice using 
the formulation of King-Smith and VanderbiltB Our cal- 
culations are both closer to the experimental values, and 
more strongly correlated with them, than the previous 
reults. However, they are still significantly lower than ex- 
periment, and the correlation is still not very good. We 
also tried including off-diagonal position matrix elements, 
and considered a modified universal tight-binding model 
that was proposed |±p incorporate non-orthogonality of 
the basis functions. Ej Unfortunately, both modifications 
were found to worsen the results. 

The details of our theoretical approach are as follows. 
We consider each zincblend compound at its experimen- 
tal lattice constant, with and without displacements of 
one sublattice along the z direction by iO.OOOlA. The 
Bloch functions are computed in the tight-binding rep- 
resentation using standard direct matrix diagonalization 
on a mesh of fc-points. According to the theory of Ref. ||, 
the electronic contribution to the polarization takes the 
form 
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where the sum runs over occupied bands and the u„k are 
the periodic parts of the Bloch wavefuctions. 
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TABLE I. The transverse charge for zincblende semi- 
conductors calculated at the experimental lattice spacing d, 
compared with perturbation estimates of Kitamura and Har- 
rison (KH) and experimental values. O- and NO- indicate 
orthogonal and non-orthogonal tight-binding models respec- 
tively, while -O and -OD refer to diagonal and off-diagonal 
representations of the position operator. 





d(A) 


O-D 


O-OD 


NO-D 




Expt.'' 


SiC 


1.88 


1.97 


2.20 


1.84 




2.57 


BN 


1.57 


1.24 


0.96 


1.01 




2.47 


BP 


1.97 


-0.09 


-0.18 


-0.23 






BAs 


2.07 


-0.39 


-0.42 


-0.54 






AlP 


2.36 


1.92 


1.61 


1.64 




2.28 


AlAs 


2.43 


1.75 


1.50 


1.50 




2.30 


AlSb 


2.66 


1.48 


1.22 


1.32 




1.93 


GaP 


2.36 


1.88 


1.57 


1.62 


0.89 


2.04 


GaAs 


2.45 


1.73 


1.47 


1.51 


0.71 


2.16 
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InP 


2.54 


2.26 


1.94 


1.99 


1.26 


2.55 


InAs 


2.61 


2.11 


1.85 


1.86 


1.07 


2.53 


InSb 


2.81 


1.86 


2.14 




0.75 


2.42 


BeS 


2.10 


1.61 


1.08 


0.71 






BeSe 


2.20 


1.56 


1.04 


0.71 






BeTe 


2.40 


1.51 


0.96 


0.57 






ZnS 


2.34 


1.89 


1.46 


0.53 


1.25 


2.15 


ZnSe 


2.45 


1.86 


1.47 


0.50 


1.15 


2.03 


ZnTe 


2.64 


2.05 


2.50 




0.98 


2.00 


CdS 


2.53 


1.98 


1.61 


1.10 




2.77 


CdTe 


2.81 


1.92 


1.53 


0.41 


1.24 


2.35 
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We are only interested in the z components of P for the 
distortions considered. After discretization in fc-space, 
these are givencl as 

where kx = {kx,ky) is discretized on a mesh of spacing 
Afc, and the contribution fcpm a string of J fc^-points 
takes the Berry-phase formM 

j-i 

(/)(k_L) = Imln ]^ det(um,kj_,fcju«,kj_,fc,+i) • (5) 

Here the argument of the determinant is a 4 x 4 matrix 
corresponding to the fact that m and n run over the four 
occupied bands. We typically use a discretization onto 
a 16x16 mesh in kj^ space, and extrapolate to J = oo 
using strings of J = 32 and J = 64 ^2 points. The trivial 
ionic contribution to Pz is added, and the value of 
deduced by simple finite differences. 

Strictly speaking, the polarization P and effective 
charge are not well-defined until the matrix elements 




Experimental effective charge 

FIG. 1. Comparison of theoretical values from this work 
(filled symbols) and from perturbation estimates of KH (open 
symbols, from Ref. 6), plotted against experimental values. 
Filled squares indicate results for compounds not considered 
by KH. 

of the position operator are specified in the tight-binding 
basis. In the context of the above formulation, these po- 
sition matrix elements are needed for the conversion (|^) 
between the M„k and ipnU- The simplest ansatz is to as- 
sume that the position operator is diagonal in the tight- 
binding representation, with elements reflecting the coor- 
dinates of the atoms. However, such an ansatz is rather 
unphysical; it would imply that the center of charge of 
an sp hybrid on an atom would lie exactly at the center 
of that atom, whereas in reality it would be displaced 
toward the principal lobe of the hybrid. We report our 
results first for the simple "diagonal" ansatz. Later, we 
discuss the effects of trying to improve upon this ansatz, 
as well as the effect of including the non-orthogonality in 
the model of Ref. |l^. |-| 

The results for the orthogonal Harrison modelQ using 
the diagonal representation of r are given in the column 
labeled "O-D" (orthogonal, diagonal) for a variety of 
zincblende structures in Table I.E3 The last two columns 
give the values of thejDerturbation estimates of Kitamura 
and Harrison (KH),EI and the experimental values, for 
comparison. The data is also represented graphically in 
Fig. 1. The filled symbols are our results; tke open ones 
are those given in Kitamura and HarrisonB (The filled 
squares represent compounds not studied in Ref. ^.) Our 
calculations shows a clear improvement, although we still 
systematically underestimate the experimental values of 
e^. The correlation between our calculations and exper- 
iment is not very good, although it should be noted that 
the lowest filled point is BN, a first row compound for 
which the model is less accurate. 

While the present results are certainly an improvement 
over the perturbation estimates of KH, there is clearly 
room for improvement. We thus investigated two possible 
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Experimental effective charge 

FIG. 2. Comparison of theoretical values using differ- 
ent variants of the tight-binding model, plotted against ex- 
perimental values. The notation is the same as in Table I. 



modifications of the tight-binding model to see whether 
they would bring the theoretical results into better agree- 
ment with experiment. First, we tried going beyond the 
artificial diagonal ansatz for the tight-binding represen- 
tation of the position operator by including some off- 
diagonal terms. Specifically, we included on-site matrix 
elements between s and p orbitals, e.g., (s|z|pz). The 
values of these matrix elements were obtained from sepa- 
rate LDA calculations on free (neutral, spin-unpolarized) 
atoms. By symmetry, off-diagonal p — p matrix elements 
of r are zero, and we assumed all off-diagonal intersite 
elements to be zero as well. The contribution of these 
extra off-diagonal terms to the polarization P was cal- 
culated as a simple expectation value, using the already- 
calculated wave functions (the Berry-phase approach is 
not needed). The results are given in the column la- 
beled "0-OD" (orthogonal, off-diagonal) in Table I, and 
are compared with the previous results (open vs. closed 
circles) in Fig. 2. Unfortunately, the correction appears 
to be in the wrong direction, and there is no apparent 
improvement in the correlation between theoretical and 
experimental values. 

Second, we attempted to improve the results by using 
a tight-binding model that includes non-orthogonality of 
the pbasis, as proposed by van Schilfgaarde and Harri- 
son.ll3 They used extended Hiickel theory to derive the 
overlap elements 



Sll'm — 



2Vw 



ll'm 



(6) 



from the original model, e; and are the onsite ener- 
gies from the same model, and K is a, parameter de- 
pending on row of the Periodic Table, chosen to fit the 
equilibrium spacings of the IV-IV crystals. The Hamil- 
tonian parameters were also renormalized following Eq. 



(11) of Ref. |T^. Some care is required in the appli- 
cation of the theory of Ref. || to this case: the inner 
product appearing in Eq. (H) has to be generalized to 
take a form like (^m.k^.fc, l^^^ j|0n,ki,fc,+i), where |(/)„k) 
is the vector of tight-binding coefficients corresponding 
to |u„k), 'S'k is the overlap matrix at wavevector k, and 
k = (kj + fcj_|-i)/2. The results are shown in the column 
labeled "NO-D" (non-orthogonal, diagonal) in Table I, 
and as the filled triangles in Fig. 2. Once again, this "cor- 
rection" is seen to act in the wrong direction, worsening 
the agreement with experiment. 

The failures of the above two attempts to improve the 
tight-binding model are disappointing, but perhaps in 
hindsight they are not surprising. For the case of the 
non-orthogonal model, a partial explanation may lie in 
the fact that the non-orthogonality was added in large 
part to improve the fit for structures that were not four- 
fold coordinated, which is not relevant here. But more 
fundamentally, we note that the model Hamiltonians we 
tested were developed by fitting to energy bands; thus, 
the fit included only information about energy eigenval- 
ues, and not the wavefunctions per se. However, the 
electric polarization is a quantity which depends sensi- 
tively on the electronic wavefunctions themselves. Thus, 
a real improvement in the tight-binding model can prob- 
ably best be accomplished by including quantities that 
are sensitive to the wavefunctions, such as e'^ values, in 
the fitting procedure itself. 

In summary, we have carried out essentially exact cal- 
culations of the transverse effective charge in com- 
pound semiconductions within Harrison's universal tight- 
binding scheme. We find a significantly improved agree- 
ment with experiment, compared with previous pertur- 
bation estimates. However, the theoretical results still 
show a systematic underestimate relative to experiment, 
by an average of 20%. Attempts to improve the agree- 
ment by including off-diagonal position matrix elements, 
or non-orthogonality of the basis, were actually found to 
lead to a worsening agreement with experiment. Based 
on this experience, we suggest that it might be helpful to 
use the effective charge as a fitting parameter in future 
tight-binding models. Such an approach might lead to 
a more accurate description of the electronic properties 
of semiconductors within this class of simple, but very 
useful, models. 

This work was supported by ONR Grant NOOO 14-91- J- 
1184. J.B. acknowledges support of ONR Grant N00014- 
93-1-1097. 
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